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Abstract 

This paper is devoted to discuss some of the features of self-similar 
(y-j ■ solutions of the first kind. We consider the cylindrically symmetric 

. solutions with different homotheties. We are interested in evaluating 

(j \ the quantities acceleration, rotation, expansion, shear, shear invariant 

and expansion rate. These kinematical quantities are discussed both 
in co-moving as well as in non-co-moving coordinates (only in radial 
direction). Finally, we would discuss the singularity feature of these 
. ^ | solutions. It is expected that these properties would help in exploring 

^ ■ some interesting features of the self-similar solutions. 
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1 Introduction 

Einstein field equations (EFEs) given by 

Rab ~ l;9abR = K>T ab , (a, 6 = 0,1,2,3) (1) 
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where R ab , R, T ab , k are the Ricci tensor, Ricci scalar, matter tensor and the 
gravitational constant respectively, are highly non-linear partial differential 
equations. To simplify these equations, we frequently impose some symmetry 
on the concerned system. Self-similarity is very helpful in simplifying the field 
equations. In Newtonian gravity or General Relativity (GR), there does not 
exist any characteristic scale. A set of field equations remains invariant under 
a scale transformation for an appropriate matter field. This indicates that 
there exist scale invariant solutions to the EFEs. These solutions are known 
as self-similar solutions. 

We study self-similar solutions because of its two important features. The 
first special feature of self-similar solutions is that, by a suitable coordinate 
transformations, the number of independent variables can be reduced by 
one and hence reduces the field equations. This variable is a dimensionless 
combination of the independent variables, namely the space coordinates and 
the time. In other words, self-similarity refers to an invariance which simply 
allows the reduction of a system of partial differential equations to ordinary 
differential equations. Secondly, self-similar solutions play an important role 
in describing the asymptotic properties of more general models. For example, 
the expansion of universe from big bang and the critical and gravitational 
collapse of a star might have self-similarity in some sense as we expect that 
the initial conditions have been lost. 

In order to obtain realistic solutions of gravitational collapse leading to 
star formation, self-similar solutions have been investigated by many authors 
in Newtonian gravity [1]. However, in GR, these solutions were first studied 
by Cahill and Taub [2]. They studied these solutions in the cosmological 
context and under the assumption of spherically symmetric distribution of 
a self-gravitating perfect fluid. In GR, self-similarity is defined by the exis- 
tence of a homothetic vector (HV) field. Such similarity is called the first 
kind (or homothety or continuous self-similarity). There exists a natural 
generalization of homothety called kinematic self-similarity, which is defined 
by the existence of a kinematic self-similar (KSS) vector field. A KSS vector 
£ satisfies the following conditions 



where h ab = g a b + u a u b is the projection tensor, a and 5 are constants. The 
similarity transformation is characterized by the scale independent ratio, a/S, 
which is known as the similarity index and can be classified into three kinds. 



£^h ab — 28h ab , 



(2) 
(3) 
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The idea of self-similarity by Cahil and Taub corresponded to Newtonian 
self-similarity of the homothetic class. Carter and Henriksen [3,4] defined 
the other kinds of self-similarity namely second, zeroth and infinite kind. In 
the context of kinematic self-similarity, homothety is considered as the first 
kind. Several authors have explored KSS perfect fluid solutions. The only 
barotropic equation of state which is compatible with self-similarity of the 
first kind is p = kp. 

Carr [5] has classified the self-similar perfect fluid solutions of the first 
kind for the dust case (k = 0). The case < k < 1 has been studied by 
Carr and Coley [6]. Coley [7] has shown that the FRW solution is the only 
spherically symmetric homothetic perfect fluid solution in the parallel case. 
Mcintosh [8] has discussed that a stiff fluid (k = 1) is the only compatible 
perfect fluid with the homothety in the orthogonal case. Benoit and Coley [9] 
have studied analytic spherically symmetric solutions of the EFEs coupled 
with a perfect fluid and admitting a KSS vector of the first, second and zeroth 
kinds. 

Carr et al. [10] have considered the KSS associated with the critical be- 
havior observed in the gravitational collapse of spherically symmetric perfect 
fluid with equation of state p = kp. Carr et al. [11], further, investigated 
solution space of self-similar spherically symmetric perfect fluid models and 
physical aspects of these solutions. They combine the state space descrip- 
tion of the homothetic approach with the use of the physically interesting 
quantities arising in the co-moving approach. Coley and Goliath [12] have 
investigated self-similar spherically symmetric cosmological models with a 
perfect fluid and a scalar field with an exponential potential. 

Recently, Maeda et al. [13,14] investigated the KSS vector of the second 
kind in the tilted case. In their recent paper [15], the same authors discussed 
the classification of the spherically symmetric KSS perfect fluid and dust 
solutions. The existence of self-similar solutions of the first kind is related 
to conservation laws and to the invariance of the problem with respect to 
the group of similarity transformations of quantities with independent di- 
mensions. This can be characterized in GR by the existence of a homothetic 
vector. 

Qadir et al. [16] have classified cylindrically symmetric static manifold 
according to their homotheties and metrics. In each case the homothety 
vector field and the corresponding metrics are obtained explicitly by solving 
the homothety equations. They found 10 solutions admitting 4, 5, 7 or 11 
homotheties. There is only one metric admitting 7 homotheties and three 



3 



metrics correspond to 4, three correspond to 5 and three correspond to 11 
homotheties. 

In a recent paper, Sharif and Sehar [17] have explored some physical 
properties of the spherically symmetric self-similar solutions of the first kind. 
This has provided some interesting features of self-similar solutions. In this 
paper, we are extending the same analysis for the cylindrically symmetric 
self-similar solutions of the first kind. The paper can be outlined as follows. 
In section 2, we shall write down the self-similar solutions of cylindrically 
symmetric spacetimes. Section 3 is devoted for the discussion of the physical 
properties of these self-similar solutions of the first kind both in co-moving 
and non-co-moving coordinates. In section 4, we shall explore the singularity 
feature of these solutions. Finally, we shall summarize and discuss all the 
results. 

2 Cylindrically Symmetric self-similar solu- 
tion of first kind 

The general cylindrically symmetric static spacetime is given by the line 
element [18] 

ds 2 = e^dt 2 - dr 2 - a 2 e x ^d9 2 - e^dz 2 , (4) 

where u, A and [i are arbitrary functions of r and a has units of length. 

Qadir et al. [17] have solved the homothetic equations and found self- 
similar solutions of the first kind (i.e. homothetic solutions). There are four 
classes of such solutions admitting 4, 5, 7 or 11 homotheties respectively. 
The class A admitting 4 homotheties has three metrics. 

The first metric is of Petrov type J, Segre type [1,111] [19] and is given 

by 

ds 2 = e 2aH ^dt 2 - dr 2 - a 2 e 2h ^d6 2 - dz 2 , (a ^ 0, 1). (5) 
The second metric is of Petrov type D, Segre type [1,(11)1] and it follows 

ds 2 = e 2ah ^dt 2 - dr 2 - a 2 e m ^d6 2 - e 2 ^dz 2 , 

(a ^(3, a, (3^0,1). (6) 

The third metric in this class is of Petrov type /, Segre type [1,(11)1] and 
takes the form 

ds 2 = e 2aln ^ } dt 2 - dr 2 - a 2 d6 2 - e 2ln ^ } dz 2 . (7) 
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The class B of homothetic solutions corresponding to 5 homotheties also has 
three metrics. The first and the second solutions are of Petrov type D, Segre 
type [1,111] and are given by 

ds 2 = e 2H ^ ] dt 2 - dr 2 - a 2 e 2h ^d6 2 - dz 2 , (8) 
ds 2 = e 2H ^dt 2 - dr 2 - a 2 d9 2 - e 2H ^dz 2 . (9) 

The third metric is of Petrov type D, Segre type [(1,1) (11)] and has the form 

ds * = dt 2 - dr 2 - a 2 e 2l ^d6 2 - e 2hl ^dz 2 . (10) 

The class C admitting 7 homotheties corresponds to only one metric which 
is of Petrov type D, Segre type [1,(11)1] and is given by 

ds 2 = e 2aH ^dt 2 - dr 2 - a 2 e 2aH ^d6 2 - e 2aln ^ } dz 2 . (11) 

This metric represents a tachyonic fluid and could be re-interpreted as an 
anisotropic fluid with an appropriately chosen cosmological constant. 

Finally, in the class, D there are three metrics admitting 11 homotheties 
other than Minkowski space. These are given by the following three metrics 

ds 2 = dt 2 - dr 2 - a 2 d6 2 - e 2H ^dz 2 , (12) 
ds 2 = dt 2 - dr 2 - a 2 d9 2 - dz 2 , (13) 
ds 2 = e 2H ^dt 2 - dr 2 - a 2 d6 2 - dz 2 . (14) 



3 Kinematics of the velocity field 

In this section, we shall discuss some of the kinematical properties of the 
self-similar solutions of the first kind, given by Eqs.(5)-(14), both in co- 
moving and non-co-moving coordinates. The kinematical properties [18] to 
be discussed can be listed as follows. The acceleration can be defined as 

U a = U a . b U b . (15) 

The rotation is given by 

u ah = u [a . b] +u [a u b] . (16) 
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The volume behavior of the fluid can be determined by the expansion scalar 
defined by 

= u%. (17) 

The shear tensor, which provides the distortion arising in the fluid flow leav- 
ing the volume invariant, can be found by 

0"a6 = «(o;6) + «(a«6) - ^Qh ab . (18) 

Shear scalar, which gives the measure of anisotropy, is defined as 

o = a ab a ab . (19) 

The rate of change of expansion with respect to proper time is given by 
Raychaudhuri's equation [19] 

^ = -\Q 2 - a ab a ab + u ab u a u b - R ab u a u h . (20) 

GST 3 

Now we discuss these properties for the above mentioned self-similar solu- 
tions. 

3.1 Kinematic Properties in Co- moving Coordinates 

First we evaluate the kinematical properties of the self-similar solutions in co- 
moving coordinates. We discuss these properties for all the classes of metrics 
mentioned above. 

3.1.1 Class A admitting four homotheties 

There are three metrics in this class given by Eqs.(5)-(7). In co-moving 
coordinates, the five quantities, i.e., acceleration, rotation, expansion scalar, 
shear and shear invariant turn out to be same for the three metrics in this 
class. The expansion scalar is zero for this class. The non-zero component 
of acceleration is 

iii = --■ (21) 
r 

The non-vanishing rotation component takes the form 

UQ1 = 2 -e aln ^ (22) 
r 
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and the shear component is 

o"oi = -woi- (23) 
Finally, the shear invariant becomes 

a 2 

- = - V ( 24 ) 

For this class, the rate of change of expansion is different for all the three 
solutions. Using Raychaudhuri's equation, given by Eq.(20), the rate of 
expansion for these solutions turns out to be 

f - £ 

dQ a 

dr r 
dQ a 

dr r 

respectively. 



— = -(1-2/3 + 3a), (26) 

— = 3 (l + 3a) (27) 



3.1.2 Class B admitting five homotheties 

For this class, all the kinematic quantities turn out to be the same for the two 
metrics given by Eqs.(8)-(9). For these two solutions, the expansion scalar 
vanishes. The only component of acceleration is 

«i = -- (28) 
r 

while the rotation component is given by 

co 01 = -e H ^\ (29) 
r 

The non-zero component of shear takes the form 

<7 i = -^oi (30) 

and the shear invariant becomes 

* = -;!■ ( 31 ) 

The rate of expansion turns out to be 

All the quantities vanish for the third solution given by Eq.(10). 
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3.1.3 Class C admitting seven homotheties 

There is only one metric in this class given by Eq.(ll) The expansion scalar 
is zero for this spacetime. The component of acceleration is 

ui = --. (33) 
r 

The only component of rotation is given by 

UJoi = 2-e ahl{ ^ ) . (34) 
r 

The shear component turns out to be 

cr i = -u i (35) 
while the shear invariant takes the form 

V -,.2 



a = -4- = Au{. (36) 



The expansion rate becomes 



dO a , 

^ = p(«+l). (37) 



3.1.4 Class D admitting eleven homotheties 

All the kinematical properties turn out to be zero for the first two solutions 
given by Eqs.(12) and (13). For the third solution, all the quantities are 
exactly similar to the quantities given for the solutions (8) and (9) except 
the rate of expansion. The rate of expansion is given by 



dQ I 

dr 



(38) 



3.2 Kinematic Properties in Non-Co-moving Coordi- 
nates 

In this section we shall discuss the kinematical properties of the self-similar 
solution in non-co-moving coordinates only in the radial direction. 
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3.2.1 Class A admitting four homotheties 

For this class, the acceleration and rotation for all the solutions turn out to 
be the same. The two components of acceleration are 

a aln(^) . Oi 

u = —e w o', m = . (39) 

r r 

The only rotation component is 

uj i = — e w o'. (40) 
r 

The expansion for the first solution turns out to be 

e = -<I±H>. (4i) 

r 

The components of shear become 

a 2aln( — ) 2 m i\ 2 7 + a 2hl(L) 

a o = 4-e V V, (Tn = — (2a- 1), o"22 = -a — — e { ->, 
r 3r 3r 

U + a ) „ (l + 7a) ^ 21n(n , , 
033 = z , O01 = e (42) 

The measure of anisotropy is given by 

1 

!> 2 

The rate of expansion turns out to be 



a = —(113a 2 - 14a + 53). (43) 



^ = _ -L(29a 2 - 2a + 14) - 2-. (44) 
rfr 9r 2V J r K J 

For the second solution, given by Eq.(6), the non-zero expansion becomes 

e = -Q£±4. (45) 

r 

The shear components are 

(T 00 = 4-e (Tn = — (a -/?), ff 2 2 = 5 <re p V V, 

r 3r 3r 

(a + 8/9) ^ln(i) 2 aln(^) /, c x 

^33 = 5 e (T i = - — (4a - j3)e (46) 

3r 3r 
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The shear invariant turns out to be 



a = — (49a 2 + 16a/? + 70/? 2 ). (47) 
The rate of change of expansion becomes 

^ = — ^(101a 2 + 9ar + 62ap + Uf3 2 + 18/?). (48) 

GST 9H 

For the third metric, given by Eq.(7), the expansion becomes the same as for 
the first solution. The non-vanishing components of shear are 

a 2aln(-^) 2 1 2 

a 00 = 4-e Kr o\ (Tn = — (2a-l), a 22 = - — (a + l)a , 
r 3r 3r 

<7 33 = -^-(a + 7) e 21n{ ^, aoi = -^(8a - l)e aIn( ^. (49) 
The measure of anisotropy is given by 

a = _L(98a 2 + 16a + 53) (50) 
and the expansion rate turns out to be 

^ = ^(l 9a 2_ 2a _ 10) _«_ (51) 

dr 9r 2V ! r K J 

3.2.2 Class B admitting five homotheties 

Now we discuss the kinematical properties for the solutions given by Eq.(8)- 
(10) of the class B. The expansion scalar remains the same for all the solu- 
tions in this class and is given by 

= --. (52) 
r 

The acceleration and rotation components are same for the metrics given by 
Eqs.(8) and (9). The non-zero components of acceleration are 

u = -e Kr o ' ui = — . (53) 
r r 

The rotation component turns out to be 

uj 01 = -e H n\ (54) 
r 
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The non-zero components of shear for the metric, given by Eq.(8), are 

4 21n(-^) 2 8 2 2ln(L) 

r dr 6r 

2 ^ ln ( — ) / r r \ 

°"33 = o-oi = ~^-e V V. (55) 

The components of shear for the metric, given by Eq.(9), have the form 

4 21n(^) 2 2 2 

o"oo = -e an = — , a 22 = -—a , 

r 3r 3r 

8 21n(-^) 7 Wi) 

^33 = -^-e <T i = -— e v n>^. (56) 

The shear invariant is same for the metrics, given by Eqs.(8) and (9), which 
turns out to be 

* = " (57) 
The rate of expansion also remains the same for the two solutions given by 
Eqs.(8) and (9) and is the following 

^ = - -. (58) 

dr 9r 2 r 

For the third solution, given by Eq.(lO), the acceleration and rotation com- 
ponents are zero. The non-zero components of shear are 

<7ll = ~, *22 = -faVO, ( X 33 = -fe 21 <), (70! = |-. (59) 

The measure of anisotropy is 



and the expansion rate for this metric turns out to be 

de 152 



dr 9r 2 



3.2.3 Class C admitting seven homotheties 



(61) 



The solution for this class has two components of acceleration. These turn 
out to be 

Ot aln(^) . OL 

u = —e w o, m = . (62) 

r r 
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The expansion is given by 



a 



= -3- = 3«i (63) 
r 



and the rotation component is 

^oi = — e r « ; . (64) 



r 



The components of shear are 



A Ot 2aln(-^) Q « 2 2aln(-^) 1 



(7 i = -2^oi- (65) 



The shear invariant becomes 



a = 30^- = 30uj (66) 



and the rate of expansion turns out to be 

f» _>„ + 2)-2 (67) 

3.2.4 Class D admitting eleven homotheties 

For the first metric in this class, given by Eq.(12), the components of accel- 
eration and rotation turn out to be zero. The expansion becomes 

= --. (68) 
r 

The components of shear take the form 

2 ° 2 ^ 21n(^-) 1 / cnN 

<7ll = —Z-, ^22 = - — , ^33 = —ZS V V, tToi = — . (69) 

3r 3r 3r 3r 

The measure of anisotropy is given by 

o- = ^ (70) 

while the expansion rate turns out to be 

dQ 56 
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(71) 



For the second solution, given by Eq.(13), all the quantities vanish which is 
exactly similar to the co- moving case. 

For the third metric, given by Eq.(14), two non-zero components of ac- 
celeration turn out to be 

u = -e ln ^\ u 1 = - 1 -. (72) 
r r 

The expansion becomes 

= -- = mi (73) 
r 



and the rotation component is 



cu 01 = -e H ^\ (74) 
r 



The shear components turn out to be 



4 21n(^) 4 1 2 

r 6r 6r 

1 ^ ln ( — ) 



The measure of anisotropy is given by 



98 



and the rate of expansion yields 

ciG 101 1 



(77) 



4 Singularities 

In this section, we shall explore the singularities of the self-similar solutions 
given by Eqs.(5)-(14). We shall use the Kretschmann scalar to find the 
singularities of these solutions. The Kretschmann scalar is defined by 

K = R abcd R abcd , (78) 

where R a bcd is the Riemann tensor. For the solution, given by Eq.(5), it 
reduces to 

2 

K = 2^-(a 2 -2a + 2). (79) 
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Since K diverges at r = hence this solution is singular at r = 0. 

For the second solution, given by Eq.(6), the Kretschmann scalar becomes 

K = — A {2a 2 fi 2 + 3[3 4 + a 4 + a 2 {\ - 2a) + 2/? 2 (l - 2/?)). (80) 
r 4 

This also shows that the solution is singular at r = 0. 

For the third solution, given by Eq.(7), the Kretschmann scalar turns out 
to be 

K = — A a 2 [a 2 - 2a + 2) (81) 
r 4 

which again gives the singularity at r = 0. 

For all the solutions of the class B, given by Eqs.(8)-(10), the Kretschmann 
scalar remains the same and is given by 

« = ^ (82) 

It is clear that K diverges at the point r = 0. Thus the solutions are singular 
at r = 0. 

For the solution in class C, given by Eq.(ll), the Kretschmann scalar 
reduces to 

a 2 

K = 48— (1 + 2a 3 -2a). (83) 
r 4 

It is obvious from here that K diverges at r = and consequently the solution 
is singular at r = 0. 

For the last class, the Kretschmann scalar turns out to be zero for all the 
solutions. 

5 Conclusion 

Self-similar solutions in GR are very important and it is interesting to discuss 
their physical features. Keeping this point in mind, we have explored some 
kinematic properties and the singularity feature of such solutions represent- 
ing cylindrically symmetric spacetime. We have discussed their properties 
both in co-moving as well as in non-co-moving coordinates (only in radial 
direction). We have explored acceleration, expansion, rotation, shear, rate 
of change of expansion and finally the singularity. 

We obtain zero expansion in co-moving coordinates for all the solutions 
while in non-co-moving coordinates we have expansion negative/positive or 
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zero depending upon the value of a and (3. There is an exceptional solution 
in the last class, given by Eq.(13), for which all the kinematical properties 
become zero in co-moving coordinates as well as in non-co-moving coordi- 
nates. In co-moving coordinates, it follows only one component of acceler- 
ation whereas there exist two components of acceleration in non-co-moving 
coordinates. The rotation component remains the same in both coordinates 
except the factor of half in non-co-moving coordinates. We see that only one 
component of shear exists in co-moving coordinates which is equal to neg- 
ative of rotation in co-moving coordinates while four or five components of 
shear exist in non-co-moving coordinates. The measure of anisotropy is given 
by shear scalar and this scalar is much larger in non-co-moving coordinates 
than in co-moving coordinates. The rate of change of expansion is zero for 
the solutions given by Eqs.(10)-(13). Thus we can say that these solutions 
are in the state of equilibrium in co-moving coordinates. All the non-zero 
properties of each solution become infinite at r = 0. 

We find that all the self-similar solutions of the first kind in classes 
A, B, C are singular at r = while in class D, the Kretschmann scalar 
becomes zero. We conclude that when we use non-co-moving coordinates, 
the quantities may become more complicated. However, we get simpler re- 
sults in co-moving coordinates. 
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